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Abstract 



We study the classical vibration problem of a chain with spring constants 
which are modulated in a quasiperiodic manner, i. e., a model in which the 
elastic energy is J2j ^ji'^j-i ~ ""i)^) where kj = I + Acos[2'Ka{j — 1/2) + 6] 
and a is an irrational number. For A < 1, it is shown analytically that the 
spectrum is absolutely continuous, i.e., all the eigen modes are extended. For 
A = 1, numerical scaling analysis shows that the spectrum is purely singular 
continuous, i.e., all the modes are critical. 
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I. INTRODUCTION 



For some years, the electronic properties of quasiperiodic systems have been extensively 
studied. In particular, a number of studies are devoted to a one-dimensional quasiperiodic 
tight-binding model [0], 

- ipn+i - ipn-i + \V{na)iJn = Eip^, (1.1) 

where V{x) is a periodic function with period 1 ( V{x + 1) = V{x) ) and a is an irrational 
number. 

The properties of the energy spectra and the eigen states of ( [1 . 1|) are generically as follows 
PJ^: When the strength of the potential A is small, all the states are extended, i.e., the 
energy spectrum is purely absolutely continuous. As A is increased, localized states appear, 
which are separated from extended states by mobility edges in the spectrum. That is, the 
spectrum has absolutely continuous parts and pure point parts, and they are separated by 
mobility edges. When A is large enough, all the states are localized, i.e., the spectrum is 
pure dense point. 

The well-known exceptions are the Harper model and the Fibonacci model 0,^. In 
the Harper model {V{x) = cos{27cx)), all the states are extended for A < 2, while all the 
states are localized for A > 2. At A = 2, all the states are critical and the energy spectrum 
is purely singular continuous. These properties are caused by the existence of the duality 
0]. It is self-dual at the critical point A = 2. 

In the Fibonacci model V^'s take two values 1 and —1 arranged by the Fibonacci sequence. 
So V{x) is piecewise constant. The spectrum is purely singular continuous and all the states 
are critical irrespective of A. 

In this paper, we study the dynamics of a one-dimensional array of equal mass particles 
connected by springs modulated in a quasiperiodic manner. The model called the modulated 
spring (MS) model is introduced in Sec. 2. A quasiperiodic case was first studied by De Lange 
and Janssen [0 . Their numerical results show the intricate Cantor-set like structures of the 
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spectra. Janssen and Kohmoto applied the multifractal method [§,0 and conjectured 



that the model has a pure singular continuous spectrum at A = 1 and has a mixed one for 
A < 1. However, their system sizes are not large enough. In fact, their claim for A < 1 is 
not correct as shown below. 

Incidentally the MS model is related to a two-dimensional tight-binding model in a 
magnetic field with next-nearest-neighbor hoppings which is introduced in Sec. 3. 



Using the result of Thouless |Tl[| for the Lyapunov exponent it is shown that all the eigen 



modes are extended for A < 1 in Sec. 4. For A = 1 numerical scaling analyses are required to 
determine the types of the spectrum. The results are presented in Sec. 5 and it is concluded 
conclude that all the eigen modes are critical for A = 1. This is presumably due to the fact 
that the spring constants can be infinitesimally close to zero. 

II. MODULATED SPRING (MS) MODEL 

The equations of motion of the one-dimensional MS model are 

u'^Uj = kj{uj — Uj^i) + kj^i{uj — Uj^i) (2.1) 

where m's are displacements of the atoms and fc's are the spring constants: 

kj = l + Acos[2n(T{j -l/2) + 9]. (2.2) 

Note that A cannot be greater than one since the spring constants /c's must be positive. Thus 
< A < 1. The MS model can be regarded as a one- dimensional electronic tight-binding 



model by rewriting with (|2.2| ) as 

— kjUj^i — fcj+itij+i -|- VjUj = Euj, (2.3) 

with 

Vj = kj + kj+i - 2 = 2Acos{-Ka)cos{2TTaj + 9), (2.4) 

and 

E = u;^ -2. (2.5) 
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III. TIGHT-BINDING ELECTRONS IN A MAGNETIC FIELD ON THE 
SQUARE LATTICE WITH NEXT-NEAREST-NEIGHBOR HOPPING 



Consider the tight-binding model on the square lattice with nearest-neighbor hoppings 
ta and tb, next- nearest-neighbor hopping tc, and magnetic flux per plaquette (see Fig. 1). 
In the Landau gauge {A = (0, By) ) the tight-binding equation is 

Since this equation is translation invariant in the y-direction, the Bloch theorem is applied 
and one may write 

V'm.n = e-*'^^™ . (3.2) 



By substituting (|37^), (lOp reduces to a one- dimensional tight-binding equation 



-{ta + 2tcCos[k + 27r0(m - l/2)]^„_i - {ta + 2tcCos[{k + 27i(f){m + l/2)])(^„+i 

—2tf,cos{k + 2'K<j)m)ipm = Eipm ■ (3.3) 

Since the original 2d problem has a symmetry with respect to an interchange of x and y 
^ has a duality. Substitute 

^^ = Y, /^.e*[^'"+(2^'^'"-'=)jl , (3.4) 



to ( p.3|) , then one gets 



-{tb + 2t,cos [K + 27r0(j - 1/2)] - (4 + 2t,cos [{K + 27r0(j + l/2)])/,+i 

-2taCos{K + 27^(1)3) fj = Efj (3.5) 

which has the same form as (13.3) • 
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IV. LYAPUNOV EXPONENT 



Using the duality in the last Sec. 3 Thouless [0] obtained the Lyapunov exponent 7 
inverse of the localization length) for (|3.ij|) as 



ln{th + Jtl - 4t2) - ln{ta + Jtl - Ul) for 2t, < < h 



7 = i ln{h + ^tl - Atl) - ln{2Q for < 2t, < h, (4-1) 

for ta <h < 2tc. 

From ( [4 .11 ) the Lyapunov exponent is strictly positive for 2tc < h and ta < t^. Thus 
all the eigen states of ( p.3| ) are exponentially localized and the spectrum is pure point. 
The duality ( |3.4| ) which maps ( p.3| ) to ( |3.5| ) by exchanging ta and ti, is essentially a Fourier 
transformation. So a localized state is mapped to an extended state. Therefore all the states 
are extended for if 

2tc < tb and ta < tt (4.2) 
Note that the MS model (|2]|) with and are the same if 

ta = 1, tc = A/2, tb = Acos(7ra)and = 6* + vr. (4.3) 



For A < 1, ( [4.2|) is satisfied. This means that all the eigen modes of the MS model \2. i| j is 
extended for A < 1. 

For A = 1 we have 2tc = ta and tb < ta- In the dual situation (2tc = tb,ta < tb) one 
has 7 = from ( |4.1| ) and the states are either critical or extended. Thus in the original 
situation the eigen modes are either critical or localized. In order to determine whether the 
eigen modes are critical or localized for A = 1 numerical analysis is required.lt is presented 
in the next section. 

V. NUMERICAL SCALING ANALYSIS 

We take a to be (V5-l)/2 (the inverse of the golden mean) and 6 to be 0. This irrational 
number cr is approached by rational numbers Fn-i/Fn {Fn-i/Fn —>■ a as n 00), where 



Fn is the nth Fibonacci number defined by F„ = + Fn-2 and Fq = Fi = 1. Thus in 
the numerical calculations a is approximated by which makes the system periodic 

with period F„. The spectrum consists of Fn bands. As n is increased, i.e., as Fn-i/Fn 
approaches to a, each band splits into sub-bands and the width of each band approaches to 
zero. We analyze scaling behaviors of the band widths for large n's. For the band width, 



we adopt the width of u"^ or E in (p^ ) and denote it as diu"^). 

For localized modes the band widths are expected to decrease exponentially with 
while for extended modes they should behave as 1/Fn [|l5l. Thus, if we define a scaling 
exponent a by ~ eigen modes are localized if a is and are extended if a is 

1. In addition, extended modes with a = 1/2 may exist. These correspond to the remnant 
of van Hove singularities at the edges of band clusters. This behavior was observed in the 



Harper model in the extended region [|T^ . For critical modes, which correspond to a singular 
continuous spectrum, a takes a continuous range of values. As n is increased by 2 or 3, each 
band splits into three subbands. Therefore, each point in the spectrum of the quasiperiodic 
limit (n ^ oo ) is specified by an infinite series of 1, and -1, which represents the upper, 
middle and lower sub-bands, respectively. In Fig.l F„5(c(j^) is plotted against n for several 
points in the spectrum for (a) A = 0.95 and (b) A = 1.0 ||TB[. For A = 0.95, the widths 
scale as 5{uj'^) ~ 1/-Fn except for the points identified by {C1C2 ■ ■ ■ — 1 — 1 — 1 ■ ■ ■} or 
{C1C2 •■■111 ■ ■ ■} where they scale as S{uj'^) ~ (1/F„)^. The latter scaling with a = 1/2 
are remnants of the van Hove singularities. The spectrum is absolutely continuous, i.e., the 
modes are extended at least for those in Fig. 2(a). 

For A = l,6{uj'^) scales with various powers of 1/Fn as shown in Fig. 2(b). It is an 
interesting result that the lowest edge band { — 1 — 1 — 1 — 1 — 1 ■■■} seems to scale exactly 
as S^uj"^) ~ (l/Fn)^. The scaling a = 1/3 is different from the edge scaling of absolutely 
continuous spectra a = 1/2 which represents a remnant of a van Hove singularity. Thus it 
seems that usual extended modes do not exist in the vicinity of = 0. 

In order to clarify the global properties of the spectrum, we perform the multifractal 
analysis P,[TO[|. We calculate /(a) of the whole spectrum which is defined by n{a)da ~ Fl^"\ 



where fl{a)da is the the number of bands whose scahng exponents a's he in an interval [a, 
a + da]. When we calculate /(a) using the statistical mechanics-like formalism |]10[ a smooth 
/(a) on [amin, c^max] is always obtained for a finite system and it is not possible to distinguish 
spectral types. Therefore it is crucial to take an extrapolation n ^ oo. 

In Fig. 3, (a) a^^l^ and (b) fj^}^ (calculated with a„ = Fn-i/ F^,) for A = 0.95 are 
plotted. It is clearly seen that Umax = 1 and f{amax) = 1 in the n ^ oo limit. Similar 
results are obtained for other values of A which are less than one. These numerical results 
are consistent with the analytical result in Sec. 4 that the spectrum is absolutely continuous 
for A < 1. 

In Fig. 4, f{a) which is obtained by extrapolation n oo for A = 1 is shown. This 
gives the convincing evidence that the spectrum is purely singular continuous rather than 
point like. Thus we conclude that all the eigen modes of the MS spring model is critical 
at A = l. 
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FIGURES 

1 The square lattice with nearest- neighbor and next- nearest- neighbor hoppings in the 
magnetic field B whose flux quanta per plaquette is (p. 

2 Fn5{u;) against n for several modes with (a) A — 0.95 and (b) A = 1.0. 

3 (a) amaxin) and (b) fmaxin) against 1/n for A = 0.95. 

4 f{a) for spectrum with A = 1. 
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